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Abstract 

Chiral symmetry breaking and mass generation is studied in a vectorial, confining asymptotic 
free gauge theory. Using the Schwinger-Dyson equation in improved ladder approximation, we 
calculate the fermion propagator in the whole Minkowski space. The estimate for f n and the 
dependence of physical mass on a coupling strength is provided. We focus on the extraction of 
spectral function of fermion propagator in the strongly coupled regime. Our calculations indicate 
that up to the crossover between walking and QCD-like running dynamics, the real pole of the 
propagator is not excluded and very likely it is actually developed at zero temperature theory. 
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I. INTRODUCTION 



Chiral symmetry plays a crucial role in particle physics. Our recent understanding of 
low energy E ~ ^qcd hadron physics is heavily based on this phenomenon. The observed 
pions are naturally identified with the approximatively massless Nambu-Goldstone modes 
of broken SU(2)l x SU(2)r chiral symmetry. The associated quark mass generation is 
driven by the QCD strong interaction of light quarks via the colored gauge boson of SU(3) C 
Yang-Mills theory. Although experimentally observable small current masses enter the QCD 
Lagrangian, the huge amount of constituent quark mass is generated dynamically. 

There are also interesting proposals to solve hierarchy problem of the Standard Model 
without point-like Higgs, which invoke strong interactions. They include Technicolor mod- 
els P, Topcolor assisted Technicolor models 0, Extended Technicolor Models 0,0] with 
walking effective coupling and even the extradimensional one In realistic versions 
of models, the part of the chiral group is a subgroup of standard model SU(2) W x U(l)y, 
while the breaking of chiral symmetry automatically implies spontaneous electroweak sym- 
metry breaking. Hence these models are very ambitious, as they should explain the observed 
masses of fermions, but they must give the correct masses to W^, Z as well. For a rather 
recent review of electroweak symmetry breaking in Technicolors see Q. 

In the above mentioned models the particles masses are mostly (if not fully) generated via 
loops. If the coupling constant is small enough we can say that the chiral symmetry plays 
the role of 'custodial symmetry' which prevents the particle from receiving a mass. Hence 
relatively strong coupling is needed, the dynamical fermion mass generation is therefore in- 
trinsically non-perturbative phenomenon of quantum field theory. There are some challeng- 
ing problems connected to the singularity structure of Green's functions in strong coupling 
models like QCD or Technicolors: e.g., questions about the presence and the position of 
propagator poles. In view of recent achievements in lattice calculations and in numerical 
results obtained by employing Schwinger-Dyson equations in the Euclidean space, solutions 
obtained in the the whole Minkowski space should bring new insights for these models. 

Let us mention some background and related work. The gap equations have been used 
for many years to study spontaneous chiral symmetry breaking in field theories [H, 0, 0, 
12L k| Solving the Schwinger-Dyson equation directly in Minkowski space has an advantage 



of providing the solution also for timelike momenta, which are not directlyaccessible in 
Euclidean formalism. Up to date, only few attempts or suggestions 3 IE 3, 03 to solve 
the problem of dynamical chiral symmetry breaking in the whole Minkowski space are known. 
In our presented approach, we assume and employ dispersion relation for the dynamical 
mass function E. However we should note that the resulting integral representation for the 
fermion propagator cannot be identified with the standard Lehmann spectral representation. 
Thus giving up the assumption of a Lehmann spectral representation for propagator, which 
seems to be unavoidable due to the high momentum asymptotic of the fermion correlators, 
we obtain perfectly stable solutions even far away from the chiral phase transition. We 
compare our Minkowski solutions with ones obtained independently in Euclidean space. 

The paper is organized as follows. In the next section we will introduce the model and 
formalism, in section IIHI we explain the solution of gap equation in Minkowski space, in 
section IIVI we will check the linear approximation employed. In section |V] physical scaling 
law will be provided, i.e., the dependence of the pole mass on the coupling. We discuss 
further possible progress in the Conclusion and Outlook. 
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II. TOY MODEL AND BASIC FORMALISM 



As our toy model consists of just one interacting quark, there is just one massless pion 
associated with broken chiral symmetry U(l) x U(l)/U(l). Once the chiral symmetry 
is spontaneously broken, the massless pole appears in the fermion-boson vertices and the 
associated current of broken generator induced non-zero transition between vacuum and the 
pion itself. In that case the pion decay constant 

<0|#(0)|tt(*) (2.1) 
can be estimated by the Pagels-Stokar formula [l| 



fl = -UN c J 



d'p nP 2 ) L, 2 , P 2 dZ(p 2 ) 



(27r) 4 (p 2 -S 2 (p 2 )) 2 



2 dp 2 



(2.2) 



where £ is mass function defined through the fermion propagator: 



wherein we take renormalization wave function Z(p) = 1 from now since its momentum 
dependence is unimportant in Landau gauge which we will use in our calculation. 

The basic input into the relation for the pion decay constant is the dynamical mass S, 
which we shall assume here to satisfy the dispersion relation: 

S(P 2 ) = Id. - P{UJ \ , . (2.4) 
J p 2 — lu + ie 
o 

In that case we approximate the E 2 in the denominator of S by the pole value m 2 , i.e. 

oo 



p 2 — S 2 (p 2 ) p 2 — m? + ie 

where the pole mass is determined self-consistently as 

S(m 2 )=m. (2.6) 

As we will see in the next section the determination of the function p and physical mass m 
is not difficult task in our denominator approximation. 

The dynamical mass is approximately constant function in the infrared p 2 < A, while it 
decreases like power of momentum 

£ * ~ 2 (2-7) 

for p 2 » A 2 in an asymptotic free theory characterized by the logarithmically behaving 
running coupling above the scale A. Recall the QCD scaling here: E ~ ^qcd{= 4:00MeV) ~ 
S(0)(= 350MeV) ~ 4/^. In our model the selfenergy behaviour ()2.7|) is achieved a bit faster 
as a consequence of the used Pauli- Villas regulator. When £ approaches A, then A represents 
the confinement scale 0] and confinement goes hand by hand with the spontaneous chiral 
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symmetry breaking. Unfortunately, the denominator approximation method of this paper 
are then no longer useful in this case. Nevertheless, the proposed denominator approximation 
should be rather good |l2| for walking (slowly running) theories when m? « A 2 . 

Here, a short but important digression is in order. Clearly, the propagator given by 
formula (J2.5)) exhibits a real-value pole. This gives rise the natural question of using a 
Lehmann representation (see the first two formulas in the Appendix) for this purpose. This 
has been attempted in the paper 17[ wherein it was found that a spectral representation with 
a single real pole part explicitly divided is applicable only with crude restrictions. Indeed, 
it gives reliable result only up to the scale A and for the couplings which are very close 
to the critical value a c . It seems that 1/p 2 behaviour of the TrS(p) (Tr is taken over the 
Dirac indices) is dictated too strongly by Lehmann representation, which contradicts 1/p 4 of 
TrS(p) decrement typically found in the Euclidean space studies. It was also shown in [TtJ 
that the propagator looses its assumed analytical property in this case. It is stated without 
proof here, that the ultraviolet asymptotic of the correlators strongly suggest the absence 
of Lehmann representation for propagators in models with the dynamical chiral symmetry 
breaking. Perhaps and with great care, the Lehmann representation can be used as a low 
energy approximation. 

In addition we neglect the terms with the derivatives and substitute our integral Ansatz 
()2.5|) into the Pagels-Stokar formula ()2.2|) for the pion decay constant. Doing this explicitly 
we get 

Sl = ~MN C I" **** 1 7^ ,, ...,± M *1, (2-8) 



(27r) 4 (p 2 — m 2 ) 2 (p 2 — x 2 )(p 2 — xl) 



After the momentum integration we arrive into the following formula: 



J IT 



4vr 2 



00 d Xl dx/ {xMX2) 
n Xi — x 2 



ln(m 2 /xi) ln(m 2 /£ 2 ) 



X 1 



x\ 



(2.9) 



This can be further approximated since one can easily recognize that up to the presence of 
log's the double integral factorizes and turns to be square of the pole masses m 2 exactly. For 
this estimate we approximate slowly varying log's in (|2.9|) bv a constant K = In (m 2 /x max ) 
where x max is the position of \p\ function global maximum. Anticipating the typical be- 
haviour of \p\ it 'explodes' from the threshold where it is zero to its maximum which lies not 
so far. Doing this explicitly and using ()2.6|) we immediately get: 



2/tt 



ni 



(2.10) 



Numerically K ~ 0.5 when S(0) ~ A which by about 50 percentage underestimates an 
experimentally well tested QCD result. It is not so bad estimate because of denominator 
approximation, Eq. (J2.9|) represents at least an interesting formula which relates the physical 
pole mass m, the absorptive part of dynamical mass function p and the decay constant of the 
Nambu-Goldstone boson. Finally notice that it is impossible to obtain such result within 
the use of Lehmann representation since the formula appears to be ultraviolet divergent in 
that case. 

Of course, in any case /„- can be calculated purely from the Wick rotated results on E 
(see Appendix ) and this is not clearly a manifest quantity where the knowledge of the 
timelike structure of the propagator is necessary. Various timelike form factors are the right 
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possible cases. To be able to calculate them from first principles in QCD the nonperturbative 
knowledge of S in the whole Minkowski is required. 

The vectorial interaction with gluons leads to the gap equation which in the Landau 
gauge and one skeleton loop approximation reads 



nv 2 ) = -w 2 c 2 (r) J 



lie 



a/3 + (p-q)a(p-q)p 

(p - q) 2 



G{p-q)S{q) lp . (2.11) 



where C^-R) = (N 2 — 1)/2N for fermions in fundamental representation of SU(N). To 
implement the running coupling, a very simple Ansatz for the product of gluon and vertex 
form factors has been used in [17]. It was chosen such that for a large Euclidean momenta 
such running coupling a run behaves like 



a. 



M 2 ) 



,(o) 



1 + 31 

A 2 



(2.12) 



which is continuous modification of sharp cutoff used in a recent studies |11| . 

Implementation of a run is easily achieved by adding negative term to the usual boson 
propagator, i.e. 

G(l) = -= 1 7T — = * , . (2.13) 



q 2 — m 2 + it 



A 2 + 



in Eq. (|2.11|) where we have introduced also an effective gauge boson mass m g . 

Further note, that we prefer to have only one dominant scale and we do not attempt to 
solve the model in all its possible peculiar phases. For this purpose we fix the ratio such 
that 

A 2 

-5 = 10. 2.14) 

m 2 g 

Throughout this paper we will express the dimensionfull quantities in the units of m g = 1 
and/or A if stated explicitly. 

Let us stressed at the end of this section that the scale A mimic AqCD or characteristic 
scale of an asymptotic free gauge theory and it has nothing to do with the renormaliza- 
tion/regularization procedure which does not take place here. Asymptotic freedom is nec- 
essary condition which make our calculation meaningful. In the real QCD the high energy 
modes are damped by logarithm which is the main difference here. As we will see, the 
running (or walking) of the effective coupling largely affects the asymptotic behaviour of the 
fermion mass function, it behaves like l/p 4 while the mass of chiral quarks behaves like 1/p 2 
for p 2 » A 2 QCD 



III. SOLVING GAP EQUATION FOR p 

The simple integral representation Ansatze allows us to convert the complicated equations 
with singular kernel to the real regular equation for the function p. The key point is the 
analytical momentum integration in the gap equation with self-consistent arrival to the 
desired dispersion relation, i.e., into the form of propagator which we have assumed and 
from which we have started the calculation. The only assumption we need is the uniqueness 
of the dispersion for p. 
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Applying the trivial algebraic identity on integral representation of S 



dx 



p(x) 



(q 2 — x + ie)(q 2 — m 2 + it) 



dx- 



p(x) 



x) 



1 



1 



m 2 + ie q 2 — x + ie 



(3.1) 



and abstracting from various prefactors, the integral kernel reduces to the differences of 
one loop perturbation theory integrals, each of them being regularized by Pauli-Villars. 
Consequently we can immediately use this result. 

The momentum integration is very straightforward since it follows standard tricks. E.g., 
the Feynman parameterization is a convenient way to perform it explicitly. At the end the 
Feynman parameter is converted into the integration variable of an appropriate dispersion 
relation. The integral is finite and does not require any subtractions as it must in the theory 
without explicit mass term. For convenience, the relevant formula can be found for instance 
in the paper [20] or in any standard textbook. 

Hence after substituting the integral Ansatz (|2.5|) into the gap equation, making a trace 
and integrating over the momentum it leads to the following expression for the dynamical 
mass: 



p{u) 
X (u;a,b) 



du- 



p(u) 



a" 



(47T) 



dx 



uj + ie 
m 2 — x 



X q (uj; m 2 , m 2 ) - X (uj; A 2 , m 2 ) 



X (u] m 2 x) + X (u; A 2 , x) 



^{uj - a-b) 2 - Aab) 



e(u-(y/a- Vbf 



(3.2) 
(3.3) 



where G represents for the usual Heaviside step function and we use following shorthand 
notation: 

C 2 (R)g 2 



a 



(3.4) 



Assuming of uniqueness DR the Eq. (|3.2|) represents a homogeneous integral equation to 
be solved. The Eq. (13.2)1 has a regular kernel. In general the homogeneous equation needs 
a suitable 'normalization condition'. Here this can be easily identified with the condition of 
on-shellness (|2.6)L i.e. 

p(x) 



dx^p— = m (3.5) 

' (m g +m) 2 Tfl — X 

where we have explicitly indicated the branch point as dictated by the arguments of step 
functions in (|3.2|) . Using Rel. (|3.5)L the equation to be actually solved numerically now 
reads: 



p(u) 



of 



(4tt) 
— / dx- 



m (X (u; m 2 , m 2 ) - X (uj; A 2 , m 2 )) 
(X (u; m 2 , x) - X (u; A 2 , x) 



(3.6) 



Using the method of iterations it is found that the integral equations is perfectly conver- 
gent and stable for almost any combinations of a* and m (without imposed (|3.5j) ). Hence 
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FIG. 1: Dynamically generated mass functions at timelike regime. Dashed (solid) line represents 
the solution for a* = 5(a* = 4). Thick (thin) line stands for dispersive (absorptive -p-) part 
respectively. 

the actual search of the true solution is very straightforward. For a given a* we scan suitably 
chosen interval of m in such a way that the identity ()3.5j) is achieved with required accuracy. 

For the purpose of comparison we solve Euclidean gap equation. After the Wick rotation 
and the angle integration it can be cast into the following form: 

a* r -m 2 g + A 2 + ^X^y^) - ^\(x,y,A 2 ) 

M(x)=m + - — / dyM(y) y - , 3.7 

8tcx Jo y + M 2 (y) 

where the positive variables x = p 2 E = —p 2 , y = q% = —q 2 and the dynamical mass is simply 
M(x) = S(p 2 ). A represents standard Kallen triangle function 

X(x, y, z) = x 2 + y 2 + z 2 — 2xy — 2xz — 2yz (3.8) 

First, let us discuss the timelike momentum behaviour. The resulting absorptive part 
of the dynamical mass function starts with negative infinite derivative from a zero at the 
threshold u 1 ^ 2 = m g + m reaching the first most robust peak. Then after going through two 
additional local maxima (minima) it behaves exactly like 

p(p 2 ) ~ ^ ( 3 - 9 ) 

forp 2 >> A 2 without any additional oscillation. The resulting mass function is displayed in 
Fig. njfor two values of the coupling a* = 4 and a* = 5. Most interestingly, to see a detailed 
behaviour we plot the absolute value of real (dispersive) and absorptive (imaginary) part 
of the dynamical mass in Fig. |21 We argue that the observed "resonant like behaviour" 
Green's functions is an inherent property of strong coupling -supercritical- dynamics. That 
this must be so for higher point Green's functions due to the formation of hadronic bound 
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FIG. 2: Log-log display of the absolute values of real an imaginary parts of dynamical mass for 
a* = 5. 

states and resonances is an experimental fact and indeed an almost trivial statement in 
QCD. Recall that the similar thing could happens to propagators of quarks is non-trivial 
observation. Recall that, this is obtained even without nontrivial inputs from higher order 
corrections involving mainly the dressed vertices which are connected directly with the 
physical spectra e.g. the lowest lying meson states tt, p, u>, a... We should stress at this 
place that there are only tiny numerical errors ( we use for integration number of 990 
Gaussian mesh points in the both cases presented, however using few hundreds with the 
suitable density is always justified). We mention here also that the principal value integral 
over the p has been calculated analytically, the result of which is presented in the Appendix 
for completeness. The cancellation among negative and positive contributions in the integral 
()2.4|) gives the correct -p~ 2 - ultraviolet decreasing log moderated behaviour of the dispersive 
part of the dynamical mass function. 

The dynamical mass function for the spacelike momenta is displayed in Fig. El Being 
below the branch point the integration washout the oscillating behaviour and the dynamical 
mass function is a smooth function. We can see in Fig. El that not so far from the critical 
coupling a* = 3.8 Minkowski and Euclidean results reasonably agree with each other. The 
numerical results of our studies point towards an analytical structure of the 'quark' propa- 
gator with a real pole and branch point singularity on the real axis. Close to the chiral phase 
transition, the suggested nature of this singularity is determined with good confidence. 

The impact of difference between walking and running SU(N) gauge theories has been 



recently investigated in Refs. 11]. Our results give additional information: vanishing dif- 
ference between Euclidean and Minkowski solution suggests that there are no additional 
complex singularities in the complex plane of the momenta. In Technicolor slang: in walk- 
ing theory the fermion propagator singularity is dominated by a real pole and associated 
particle mode can be unconfined. Disappearance of the real propagator pole would be a 
certain indication of confinement in any model which is heavily based on the walking be- 
haviour of the coupling constant. Recall the possibility of chiral symmetry breakdown and 
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FIG. 3: Dynamical mass function in the spacelike regima of momentum. Spectral Minkowski 
solutions are labeled by SM , the solutions obtained in the Euclidean space are labeled by EU. 




FIG. 4: Scaling laws as it is seen in the Euclidean and Minkowski approximations. 



the absence of confinement is supported by an arguments in 2 • 

In Fig. [H we show the solution for dynamical mass as a function of a*. In the walking 
regime (1 — 1.1) cv* it can be parametrized by the known 0,E3| scaling low 



S(0) ~ C ■ Ke -[K<**/<-D-^] j 



(3.10) 



where C = 0(l)and k ~ 7r . 

Increasing the coupling and leaving thus the walking coupling domain and slowly ap- 
proaching QCD limit S ~ A a certain discrepancy appears. In present, the nature of the all 
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FIG. 5: Checking denominator approximation as described in the text. The coupling strength is 
a* = 5 and the obtained pole mass is m = 0.695m 9 

amount of the discrepancy cannot be determined with confidence. However we point out the 
main part is due to the linearized approximation (|2.5J1 . The other, more speculative reason, 
could be that the pole of the propagator, located for subcritical case on the real axis, moves 
gradually into the complex plan receiving non-negligible imaginary part. Hence the integral 
representation employed for our calculations would become inappropriate. 

IV. CHECKING THE APPROXIMATION fl220) 

At present we do not know any reason why the spacelike Euclidean calculation should dif- 
fer from the 'true' Minkowski space calculation on the real negative axis of four-momentum. 
Hence the calculation performed in the Euclidean formalism can serve as a good guidance. 
To observe how much the neglect of a running mass in the denominator (J2.5j) affects the 
behaviour in spacelike domain, we make a similar approximation in the Euclidean space 
counterpart. For this purpose we replace in Eq. (|3.7|) the appropriate part of the propaga- 
tor follows 

M(y) ^ M(y) 
y + M 2 (y) y + M 2 (0) 1 ' ' 

hoping that M(0) lies sufficiently close to the true pole value m. Note, albeit we have 
estimated m in our Minkowski treatment we cannot use this value, since this is already 
underestimated just due to the discussed approximation. 
The solution is shown in Fig. El for a* = 5. 

The corresponding approximation is labeled by the dashed line and much better agree- 
ment with the Minkowski approximation is achieved. It leaves little space for other effects 
and indirectly indicates that the position of propagator singularities actually lies on the real 
axis. We should mention that the approximation ()2.5|) exhibits a similar effect for all studied 
value of a*. On the other side, the splitting of the obtained scaling law functions becomes 
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FIG. 6: Comparison of physical scaling law with S(0) in Minkowski approximations. 



rather large when m/A ~ 1. To avoid some speculations we conclude that it simply calls to 
go beyond the approximation employed. 



V. PHYSICAL SCALING LAW 



The pole mass of the propagator can be the only physical observable here. In this section 
we will largely simplify and speculate. We disregard the systematic error which follows 
from the approximations and assume the ratio m/E(0) being reliably estimated. We leave 
aside the complicated task of confinement-deconfinement-chiral phases 0, 21 1 which lies 
beyond the scope of this paper and certainly beyond the ability of the toy model used here. 
We will simply assume that the modeled propagator is a part of the theory in which the 
singularities of the Green's functions persist in the S-matrix (That the opposite can happen 
is well known). In this case the pole mass of the propagator is the only possible physical 
observable here. In a gauge theory, the value of dynamically generated m cannot be affected 
by a gauge fixing choice and cannot depend on the associated auxiliary parameters. To make 
some approximations in the infinite gap equations system is necessary, e.g., the truncation 
of the equations system, have always some unwanted impact on the calculation scheme 
(in)dependence. Many studies have investigated this important, but complicated issue. For 
some discussion in the context of Euclidean Schwinger-Dyson formalism see [9|. The value 
S(0) is usually believed be an approximately scheme invariant if its value is rather close to 
the physical m. Then there follows important message from the Minkowski study attempted 
here. We present the physical scaling low in Fig. El and compare with the more standard 
infrared value S(0). The observation of large splitting of the lines in Fig. El is clear. The 
ratio m/S(0) deviates from the save value one rather early. In fact, approximating m by 
the infrared S(0) is justified only for slowly running theories. 
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VI. CONCLUSION AND OUTLOOK 



In the ladder approximation of gap equation with the running coupling modeled by Pauli- 
Villars regulator we obtain a dynamical mass function in the whole Minkowski space. The 
momentum space gap equation has been converted to the real regular equation for a real 
function p. The procedure is for us surprise perfectly stable for all couplings studied. We 
compare the solution with the ones obtained in the Euclidean formalism. The employed 
approximation works rather well for the coupling not so far from the critical one. The 
Minkowski results appears to be underestimated by about 50 percentage when we get scaling 
£(0)/A = 1. The physical pole mass has been compared with the infrared one and the 
physical scaling has been established. 

We did not plot the the ratio f n /m as a function of a * (, A). Using a rough estimate by 
the Pagels-Stokar formula we found that up to deep walking limit this ratio has only weak 
dependence on the coupling strength a*. It decreases from the value f n /m = 0.5 when a is 
few percentage above a* to the value 0.25 when a ~ 10 (QCD limit). 

The absorptive and dispersive parts of the dynamical mass function are not positive 
definite functions in the timelike regime, while the later decreases monotonically down in 
the spacelike regime in a shape already known from the Euclidean studies. The resonant 
behaviour of Green's function in the timelike regime can be very inherent to the model with 
strong coupling, especially for the light quarks in real world QCD. If such behaviour is the 
reality and general feature of the fermion propagator in chiral breaking phase it calls for a 
better understanding of its physical impact. The presented toy model hopefully provides a 
reliable result for slowly running coupling. The application to many questions of electroweak 
symmetry breaking in Technicolor models is left to a future work. When one moves far away 
from walking behaviour, one loses a simplifying feature of our calculation, namely that the 
Pauli-Villars form factor in use does not reliably approximate QCD running coupling. A 
more elaborated approach suitable for such situations is presently in progress. 
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APPENDIX A: USEFUL FORMULAS 

Kallen-Lehmann representation: 

For fermion propagator in parity conserving theory reads 

S(P) =i>S„(/) + 1 = I du, • (AD 

Jc p z — uo + le 

hence assuming real pole, the only possibility is to assume Dirac delta distribution as the 
lowest mode. In that case we can write for instance for S s 

S S (P 2 ) = 2 " ■ + r - ^ , , (A2) 

p l — m 2 + xe J( m +m g ) 2 p l — uj + ie 
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where r is a residuum, which does not vanish when dynamical chiral symmetry breaking is 
attempted within this formalism. Recall the Ansatz (jA2|) works rather well when it is used in 
the Schwinger-Dyson equations' formalism with explicit mass terms and for the subcritical 
couplings 0, 0, Q . 



Dispersion relation integral 



There are few convenient possibilities how to perform integration in the dispersion relation 
for the dynamical mass. Most simple is to use numerical data for p and integrate numerically. 
This is the way which we actually follow below the branch point. On the other hand, one 
can avoid numerical principal value when evaluating M for timelike p 2 > (m + m g ) 2 by 
analytical integration. For this purpose we can use the equation ()3.6|) and substitute it into 
the (j2.4j) . arriving thus to the following regular integral: 



ReZtf) = ~ fdx^ 
Am J m 2 — 



x 



J(p 2 , m 2 , m 2 ) - J(p 2 , m 2 , A 2 ) - J(p 2 , x, m 2 ) + J(p 2 , x, A 2 ) 



(A3) 



where we have defined 



J(p 2 ,y, z) 



P 2 



7T p*-y-z 
— h arctq . — 

2 . 



ln(16yz) 



In 



P 



p —y — z + 



T-p 2 



p z 



with the following abbreviations: 



y — z + \ A 



+ e(A £ K^ ln 

p2 



-y-z + 



A 



-y - z + V\q 



(A4) 



Ao 



\y~ z \, 



t = (Vy + ^z) 2 . 



(A5) 



These expressions have been actually used for the numerical evaluation of M in the presented 
paper. 
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